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The thermal conductivity k of the quasi-two-dimensional (Q2D) organic superconductor A- 
(BETS^GaCU was studied in the magnetic field H applied parallel to the Q2D plane. The phase 
diagram determined from this bulk measurement shows notable dependence on the sample quality. 
In dirty samples the upper critical field H C 2 is consistent with the Pauli paramagnetic limiting, and a 
sharp change is observed in n{H) at H c2 \\. In contrast in clean samples H C 2{T) shows no saturation 
towards low temperatures and the feature in k.(H) is replaced by two slope changes reminiscent 
of second-order transitions. The peculiarity was observed below ~0.33T C and disappeared on field 
inclination to the plane when the orbital suppression of superconductivity became dominant. This 
behavior is consistent with the formation of a superconducting state with spatially modulated order 
parameter in clean samples. 



PACS numbers: 74.25.Fy, 74.25.Dw, 74.70.Kn 
I. INTRODUCTION 

In most of superconductors, the suppression of super- 
conductivity by a magnetic field H is caused by the 
increase of diamagnetic energy. Since the effect is or- 
bital in nature, the upper critical field H C 2 becomes high 
when the orbital motion is suppressed. Under this con- 
dition, the destruction of spin-singlet state of the Cooper 
pairs may become dominant mechanism of SC state sup- 
pression, with the H C 2 determined by the Pauli para- 
magnetic limiting field HpU. Both of these mechanisms 
give saturation in thelemperature T dependence of H C 2 
at low temperaturesou. In a number of recent experi- 
ments, however, a nonsaturating behavior of the H C 2(T) 
was observedB, and its possible relation to the formation 
of the inhomogeneous superconducting state was consid- 
ered. 

The existense of inhomogeneous state was predicted 
back in the sixties by Fulde and FerreH and Larkin and 
OvchinnikovQ, who pointed out that the stability of the 
superconducting (SC) phase can be increased above H p 
by pairing electrons with different momenta. In their 
model at low temperatures the transition from the nor- 
mal state at H C 2 proceeds into a new SC phase with non- 
zero momentum of the Cooper pair and spatially mod- 
ulated order parameter, abbreviated recently as FFLO 
state, and becomes second order, contrary to the first or- 
der transition expected at H p . Simultaneously at lower 
magnetic field an additional phase boundary appears 
within the SC domain, associated with the transition 
from FFLO state to usual SC state with zero momen- 
tum of the Cooper pair (hereafter we call it BCS state! 

Despite quite long history of theoretical studiesQiiZI, 
the experimental observation of FFLO state is contro- 
versial. The state can be formed only if orbital motion 
is strongly suppressed, so that respective upper critical 
field H C 2orb becomes larger than Hp£3, while the super- 



conductor remains in a clean limitu. These conditions 
are not met simultaneously in usual superconductors, in 
which the orbital effect is suppressed by alloying for ex- 
ample. On the other hand, favourable situations for the 
formation of thejiFLO state are found in heavy fermion 
superconductorsEj and in quasi-two-dimensional (Q2D) 
superconductors under a magnetic field parallel to the 
Q2D planeQ; in these cases, due to the large effective 
mass, H C 2orb becomes very large even in clean samples. 

Indeed in several experimental studies, pronounced 
anomalies that were observed close to the Jii 
Pauli paramagnetic limiting, range in heavy- fermionEj CJ, 
organicOtJ and cupratecJ superconductors were dis- 
cussed in a context of the formation of FFLO state. It 
should be noticed, however, that in all these studies one 
cannot rule out another effects appearing near H C 2, espe- 
cially in the resistivity p and magnetization M measure- 
ments. These can include anomalies related to the trans.- 
formations in thcj-siortex stateEH, such as .peak effectES, 
commensurabilityB3 and melting transitioned, and the ef- 
fects of filamentary and surface superconductivity^ on 
measured H C 2(T) line. Identification of the FFLO state 
through the|-alaape of the H C 2 line with upturn at low 
temperaturecj'ES is not reliable as well, since the-similar 
shape finds explanation in a number of scenariosEdi^ji 

A key property for the identification of FFLO state is 
its sensitivity to disorder. While the FFLO state should 
be rapidly suppressed by disorder, this is not true for 
the transformations in the vortex lattice. Indeed, al- 
loying was found to. increase anomalies close to H C 2cB, 
attributed initiallyc2l to the FFLO state in CeRu2 . How- 
ever, the effect of disorder is not, unique for FFLO state, 
since it is similarly importaniEJ for several models ex- 
plaining both high upper critical field values and upturn 
in the H C 2(T) line on going to T=0. 

In view of this it is desirable to make a more detailed 
experimental study on the subject. As useful technique 
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to identify the FFLO state, we adopted the measure- 
ment of thermal conductivity k. Its main advantages are 
bulk nature of the signal and insensitivity to the trans- 
formation in the vortex system. The former allows us to 
disregards filamentary and surface contributions, making 
dominating contributions to the resistive measuremcnto. 
The latter is related to the absence of Lorenz force act- 
ing on vortex lines, since the thermal current of quasi- 
particles (QP) does not create charge flow. As a conse- 
quence, the flux creep phenomena, determining p and M, 
are ppt essential in k, allowing determination of intrinsic 
if c2 E3. Within the mixed state, the electronic part of k 
in conventional SC is caused by QP tunnelling between 
vorticesa-cil, while in unconventional SC-the main contri- 
bution comes from the delocalized QPc3. Both of these 
mechanisms are not influenced by the transformations of 
the vortex state and are not accompanied by any feature 
in n(Hp3. Therefore, the transition from the BCS state 
to FFLO state, -udaich is expected to be accompanied by 
the release of QPliia, should be directly seen as an increase 
in k(H). It is also advantageous that k can distinguish 
first— <p*der transition, leading to a jump in k(H) near 
iJ c2 c3'c3 ) from second order transition giving an anomaly 
in the slope of k(H)u. 

In this article we report the study of thermal conduc- 
tivity at low temperatures in the Q2D organic supercon- 
ductor A-(BETS)2GaCLi in H aligned precisely parallel 
to the conducting plane so as to suppress orbital mo- 
tion. This superconductor is a non-magnetic version of 
A-(BETS)2FeCi4 showing field-induced superconductiv- 
ity due to a Jaccarino-Peter effectoti9. It was proposed 
theoretically rtiiat the FFLO state may be formed in both 
of these saltsE2lEil. This study became possible by notable 
technical advancement in thermal conductivity measure- 
ments in high and oriented magnetic fieldsu. We show 
that the dependence of H C 2 on temperature, inclination 
angle and sample quality, together with sharpening of 
transition in k(H) curve with disorder, is consistent with 
the formation of FFLO state in clean samples of this or- 
ganic superconductor. 



peratures. The same contacts were used for resistance 
and thermal conductivity measurements. The direction 
of the heat flow was along the longest dimension, corre- 
sponding to the crystallographic c axis within conducting 
plane. 

The cell with the thermal conductivity measurement 
unit was rotated in 3 He ambient by a double-axis go- 
niometer in a superconducting solenoid. The parallel 
alignment of a magnetic field to the Q2D plane of the 
sample was performed by measuring sample resistivity 
as a function of an inclination angle O in a magnetic 
field close to H c2 \\ and determining the position of the 
minimum in p(O) curve. The accuracy of this alignment 
was determined by the sharpness of the feature and was 
typically ±0.1°. 

Single crystals of A-(RETS)2GaCLi were grown by the 
electrochemical methodEj. The samples had typical size 
of 2x0.1x0.1 mm 3 . We studied 6 single crystals from 
3 different batches, their properties are summarized in 
Table fl A notable problem with measurements on these 
samples comes from resistance jumps on cooling. Despite 
use of soft Pt wire suppoijt-|and a precise control of the 
cooling rate at 0.5 K/minE3, resistance jumps occurred 
in every second sample, especially in the range of struc- 
tural transformation near 100 K, giving notable scatter 
in the residual resistivity value po at temperatures just 
above the onset of the superconducting transition. The 
frequent jumps were to some extent caused by unfriendly 
stress conditions in thermal conductivity unit. Contrary 
to the resistivity, effect of the jumps was not pronounced 
in the thermal conductivity. The jumps are usually re- 
lated to the formation of structural domain wallsEj, but 
this experiment indicates that the wall has different ef- 
fect on charge and heat flow. The thermal transport 
is not influenced much by the boundary provided that 
the acoustic impedance is not changed substantially and 
therefore the break in electronic heat flow at a wall is 
healed by the phonon contribution. 



III. RESULTS 



II. EXPERIMENTAL 

Measurements of thermal conductivity were per- 
formed using a standard steady state one-heater-two- 
thermometers technique. A sample was attached via an- 
nealed 10 Pt wires to a copper base-iof a miniature 
vacuum celE3 and to Ru0 2 chip resistors^ acting as both 
heater and thermometers. The thermometers were cal- 
ibrated in the temperature range 0.3 to 10 K at a set 
of magnetic fields and the magnetoresistance correction 
was made for in-field thermometry. The contacts to the 
sample were made by gold evaporation on fresh sample 
surface and subsequent wire attachment to a pad with 
Dotite silver paint. This technique allowed us to get con- 
tact resistance in the range of 100 to 500 mfl at low tem- 



At room temperature, the resistivity p of the samples 
studied was in the range from 140 to 190 fiflm, with 
± 10 % error bar coming mainly from the uncertainty 
in the determination of the geometrical factor on thin 
samples with parallelogram cross-section. The tempera- 
ture dependence of p above 90 K was very similar among 
the samples fmm three batches studied, matching previ- 
ous reportsEjQ. On cooling, p remained almost constant 
down to ~200 K, increased slightly to ~ 90 K and then 
decreased rapidly all the way down to the superconduct- 
ing transition. The magnitude of the decrease below 90 
K was notably dependent on sample batch, with varia- 
tion in the residual resistance ratio from 20 to 150. In 
Fig. [l] we show p(T) curves for two best samples showing 
no jumps on cooling. The samples belonged to batches 
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A and B, characterized by a notably different residual re- 
sistance. The samples from batch C were of high quality 
as well, however, they were notably thinner and we have 
not succeeded in getting jump-free resistance curve for 
sample #6. Worth of noting that p(T) dependence just 
above the transition temperature T c in the high RRR 
sample is well described by Independence ascribable to 
electron-electron (ee) scattering. 

The k(T) /T in the SC region is shown in Fig. |. It 
increases on entering the SC state, and then shows the 
maximum at temperature T max with magnitude strongly 
dependent on sample quality. TJiis peak is known in 
a wide range of superconductors! It is caused by the 
condensation of normal electrons due to formation of a 
superconducting gap. In the normal state, conduction 
electrons are_acting as scatterers for both heat carri- 
ers, electronsEj (via electron-electron (ee) scattering) and 
phonons (via phonon-electron (ep) scattering), determin- 
ing their mean free paths, l e and l gi respectively. On con- 
densation the mean free path increases until reaching a 
limit determined by impurity and boundary scatterings, 
therefore the difference in the magnitude of the increase 
provides another way for characterisation of sample qual- 
ity. The ratio of k/T at T c and at T max is shown in Table 
H The values are well reproduced between samples from 
the same batch. There was a clear correlation between 
the residual p above T c (when this quantity was not in- 
fluenced by resistance jumps) and the magnitude of k/T 
increase, as can be seen for the samples #1 and #2. 

The T cmic i was determined from the mid point giving 
1 /2 of po evaluated at the junction point of the tangential 
lines for p(T) in the normal and transition states. The 
derived values are approximately the same for all samples 
studied, (see Table ||). The resistive transition is rather 
broadEJ, but its width decreases systematically with the 
decrease in po- For sample #1 the resistive T cm id appear- 
ing at 5.2 K corresponds to the onset of k/T increase. In 
contrast, the superconducting transition for #2 did not 
give clear anomaly in k/T(T) at T c , which makes the 
determination of T c and H C 2 from k measurements im- 
precise at high T. However, k(H) gives clear anomaly at 
H C 2 at low T, caused by a rapid change of K e , as seen in 
both samples (Figs. H]4|), enabling high precision in H C 2 
determination at low T. 

In Fig. H we show k(T) at 0.3 K, the base tempera- 
ture of our experiment. The dependence was taken as 
a function of H of different orientation with respect to 
the heat flow and conducting plane in sample $1. The 
k(H) dependence in the field applied perpendicular to 
Q2D plane is usual as for most of the superconductors 
at low temperatures. On the field increase, k decreases 
gradually then takes a minimum followed by the rapid 
increase towards H C 2- In the normal state above H C 2, 
k decreases with field due to the magnetoresistance ef- 
fect on K e . The shape of the curve can be easily under- 
stood if we recall that k contains contributions of two 
heat carriers, electrons K e and phonons K g , with k = K e 
+K g . Since the density of phonons is not influenced by 



a magnetic field, the effect of magnetic field on K g comes 
solely from the scattering. The phonons are effectively 
scattered by conduction electrons, and as we pointed out 
above, their mean free path is notably increased within 
the SC state. In the field the phonons become scattered 
by vortices. This effect is responsible for initial decrease 
of k under the perpendicular field, as shown in Fig. |^. 
The decrease-of K g with H is gradual and is well described 
by 1/H la-wea. K e gradually increases with field at low 
temperatures (bo th , f or conventional and unconventional 
superconductors )r 4 H 61 l, while it increases rapidly towards 
H C 2 in both cases due to the intervortex tunncllingaH. 
Finally at H C 2 the increase of the electronic contribution 
is ceased and K e decreases due to magnetoresistance in 
the normal state. 

Assuming the validity of the Wiedemann-Franz (WF) 
law, we can determine K e from p measured with the same 
geometrical factor when the same contacts are used. This 
estimation becomes evidently invalid when the sample 
shows resistance jumps, since weak links at domain walls 
behave differently for charge and heat transport. In the 
samples which showed no resistance jumps this estima- 
tion gives K e equal to 30% of total k in sample #1 and 
10% in sample #2. We show these values with double 
ended arrows in Figs. || and |]. Another way to esti- 
mate the electronic contribution is to assume the same 
change of k and p with magnetic fieldE^I. This way of 
estimation produces the values which are quite close to 
the estimation via the WF ratio. 

Since K g is gradually decreasing with H and all in- 
crease in k in the proximity to H C 2, which is of most 
interest for us, is electronic in origin, we consider phonon 
contribution as a background. The justification for this 
can be found in the dependence of k on the field orien- 
tation with respect to the heat flow direction. While the 
phonon conductivity is slightly varying with the field di- 
rection, the increase near H C 2 does not depend on it. The 
magnitude of increase in k near H C 2 can be regarded as 
a lower bound of the electronic thermal conductivity. As 
can be seen from Fig. || the magnitude of the increase is 
in quite reasonable correspondence with the estimation 
of K e via the WF law. 

In the parallel field, k{H) shows basically the same 
features as in the perpendicular field, except for several 
special points. The decrease of K g is still gradual (as 
can be more easily seen in k{H) at higher temperatures, 
Fig. [|) . Due to the anisotropy of the coherence length, 
the cross-section of phonon scattering by vortices in the 
parallel-field is notably smaller than in the perpendicu- 
lar fielcH, and the decrease of the phonon contribution 
is smaller, also. The cross-section of scattering by vor- 
tices slightly depends on the orientation of the field with 
respect to the thermal flow, as can be seen in Fig. ||. 
However, K e shows notably different behavior in the par- 
allel field, as compared to the perpendicular field, aud 
this difference is directly related to the Pauli limitingED. 
When the transition at H C 2 is caused by the orbital effect, 
vortices gradually fill the volume of the superconductor. 
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On approaching H C 2 the distance between vortices grad- 
ually decreases, giving high probability of quasi-particle 
tunnelling between vortices. Because of this, K e shows 
gradual and rapid increase near H C 2- In case of Pauli 
limiting the vortex matter does not fill all the volume of 
the superconductor at H C 2, which is a consequence of a 
first order transition at H p . As a result, the change of 
K e becomes stepwise, provided the system is clean. Ac- 
tually the jump-like feature was observed in the Raidi 
paramagnetic limiting range in UPtaEl and CeCoIn 5 EI. 

It should be pointed out that the behavior of thermal 
conductivity in parallel field at low temperatures is well 
reproduced between different samples in the same batch 
and correlates with the magnitude of the peak in n(T) /T. 
In Fig. H we show the field dependence of K in parallel 
field in the samples under study, and in Table | we sum- 
marize the parameters of the n(H) curves at the base 
temperature. These parameters include the value of a 
field H* where n starts to increase, the magnitude of k 
increase from H* to H c ii Ak= k(H c2 \\)- k(H*), and the 
value of H c2 itself. 

Worth of special note is the transformation of this be- 
havior with disorder. In the low quality samples #2 and 
#5, k shows rather sharp increase in the proximity of 
H C 2, reminiscent of a broadened first order transition. 
In the high quality samples #1,#3,#4 and #6 the fea- 
ture at H C 2 greatly broadens and is replaced by two slope 
changes, with the increase of k starting at lower and fin- 
ishing at higher H than in the dirty samples. 

To elucidate the origin of this behavior we studied its 
transformation with temperature and orientation of mag- 
netic field both on inclination to conducting plane and on 
rotation parallel to it. Since these measurements are time 
consuming and the effect of the resistive jumps can not be 
completely ignored, we made the thorough studies on two 
best samples #1 and #2. In both cases the anomalous 
behavior near H C 2 disappears and k(H) becomes similar 
to the perpendicular field case for inclinations O above 
-5°. 

In Fig. U we show n(H) for the field applied parallel 
to the plane (0=0°) at several T, with the curve for the 
perpendicular field 0=90° at 0.3 K, as a reference. As 
already pointed out, the anomalous behavior near H c2 is 
observed only in the parallel field. (Fig. |J) 

IV. DISCUSSION 

As can be seen from Table [jj the behavior of thermal 
conductivity is consistent between the samples from the 
same batch, supporting that it is intrinsic. We will dis- 
cuss the results by sticking to the behavior of two sam- 
ples, for which we can get reliable value of p at low tem- 
peratures due to the absence of resistance jumps. For 
these samples the difference in residual resistivity po was 
the largest (Fig. |l|) with p—8xl0~ 7 f2m for sample #1 
and 12xl0~ 6 f2m for sample #2. In standard conductiv- 



ity theory this gives scattering time r of 0.16 ps (#1) and 
0.01 ps (#2). The l e can be estimated from the Fermi ve- 
locity averaged over the hole pocket of the Fermi surface, 
!i F «3* 10 5 m/s,N as 48 nm (#1) and 3.2 nm (#2). 

In Figj^ we show the H versus T phase diagram for 
#2 (coinciding within experimental accuracy to that of 
#1) determined from p(T) measurements together with 
the data determined from k(H) at low T. The slopes 
of the H c2 (T) near T c for the perpendicular and parallel 
fields are equal to 0.45 T/K and 4.3 T/K, respectively. 
This gives the values of the coherence length £j_(0) = 1.2 
nm (out-of-plane) and £||(0)=11.6 nm (in-plane). The 
comparison of £_l(0) (effective for H c2 for 0=0°) with 
the l e shows that #1 is in the clean limit, while sample 
#2 is on the border between clean and dirty limits. 

Within the SC state the resistivity is determined by 
the vortex motion and the increase in p starts below H c2 , 
as determined from n(H) in #2 (Fig. |^). However, in 
#1 the increase of k at 0.3 K starts at H giving p — 
0. This implies that the increase is not related to the 
vortex motion, but reflects the QP flow in the bulk. In 
the inset of Fig. |^ we show the H c2 data determined from 
the k(H) measurement shown in Fig. ^ for 2 samples 
together with data points of the onset of k(H) increase 
in $1. It is noteworthy that in the high quality sample 
the linear increase of H c2 on cooling correlates with the 
anomalous behavior of k(H). 

The anomaly observed in n(H) at H C 2 in the low qual- 
ity sample #2 is reminiscent of the behavior in the Ejj4i 
paramagnetic limiting range in UPt3 and CeCoInsOufj, 
in which the first order transition at H c2 leads to a jump 
in k. Indeed, a clear saturation of H C 2(T) at a value 
not far from the expectation for H p in the weak coupling 
BCS model (H p =l.MT c with H p in T units and T c in 
K units), can be seen below ~2 K (Fig. |). Taking H p 
as equal to the H c2 of #2 we can verify whether con- 
ditions for the formation of FFLO state are met m our 
experiment. With ff C 2or&=0.71^§^ T c =16.3 TO this 

leads to 0(= ) =2.010, which is suitable for the 

formation of FFLO state. 

Although there is no jump in n(H) as well as hys- 
teresis with field, which can confirm the first order type 
of transition directly, the transition is much sharper in 
dirty samples. On the contrary, in the high quality sam- 
ple, k(H) in the region of H c2 is broadened and replaced 
with two slope changes. This broadening in clean sam- 
ples shows that some new process is involved. Indeed, it 
is natural to expect broadening of the phase transitions 
with disorder. This is really observed as broadening of 
the resistive transition near T c in zero field, but is strictly 
opposite to the behavior near H C 2- To characterize this 
feature we plot in the inset of Fig. |^ the field H* , where 
the increase in k(H) starts in the high quality sample, 
as a function of T. Although the choice of this point at 
high T is somewhat ambiguous due to steep phonon back- 
ground variation, it is clear that H * (T) line extrapolates 
not far from the point where the H C 2 lines for the two 
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samples begin to deviate, as expected for FFLO state. 
On the other hand, in sample #2 the feature remains 
sharp, although due to much large contribution of K g , 
it is difficult to give any quantitative characterization of 
this fact. 

In Fig. H we show dependence of the H C 2 and of 
H* of sample #1 on the orientation of the field within 
the two-dimensional plane. As it can be seen, none of 
the lines shows-notable anisotropy, contrary to resistive 
measurements^], supporting their relation to Pauli para- 
magnetic limiting. 

Summarizing our experimental observations we can 
state the following. (1) The shape of the phase diagram 
in the parallel magnetic field, H c2 \\(T), depends on the 
mean free path of samples. In clean samples the line 
shows no saturation at low temperatures, while it shows 
clear saturation in dirty samples. (2) The increase of k 
near H c2 \\ gets essentially sharper in dirty samples indi- 
cating a tendency towards the first-order transition, while 
the second-order transition is observed in clean samples. 
(3) The increase of k(H) above phonon background starts 
at lower fields in clean samples, indicating an additional 
phase boundary within the SC phase at low temperatures 
and high fields. (4) The difference is clearly observed be- 
low ~1.7 K, i.e. ~ 0.33T C . (5) The anomalies are specific 
to the parallel field direction and disappear on the field 
inclination beyond O c ~ 5°, where the orbital motion 
becomes dominating. 

All these features are in line with theoretical predic- 
tions for the FFLO state. An important issue, however, 
is whether this explanation is unique. Several models 
were put forwardEj l3 to explain an— unusual upturn of 
H C 2(T) at low T in high-T c cupratesEa. These are, how- 
ever, not specific for the parallel field orientation. Among 
the models based on the low-dimensionality of the sys- 
tem, .the most relevant is the model by Lebed and Yamaji 
(LY)cJ. This model is a generalization for Q2D caser-of 
the original idea, developed by Lebed for Q1D systemsEa, 
and studied latee,as well by Dupuis, Montambaux and Sa 
de Melo (DMS)E3. We would refer to the LY theory, since 
its dimensionality is in accordance with the present case. 
The model takes into account quantum corrections to the 
orbital motion, becoming important if the length scales 
of the orbital confinement (interlayer distance d) and of 
the superconductivity are comparable. The FFLO 
and LY theories make similar predictions with respect to 
the suppression of upturn of the -EL2 on inclination and 
on increase of impurity scatteringcJ. There is, however, 
notable quantitative difference. An upturn of H C 2(T) 
in LY theory starts below T* w 0.1(d/£i(0))T c (0). In 
A-(BETS) 2 GaCl 4 this gives T* « 0.6 K, notably lowex 
than ~1.7 K. In the FFLO state the upturn is predictedEl 
to start at Ti—0.56T C , ~2.7 K, if the orbital motion is 
negligible, which is 1.5 times higher than extrapolated in 
our experiment. The account of the orbital effect, how- 
ever, is expected to shift the transition point to lower T, 
improving matching with the experiment in the present 
case. 



In addition, the LY, Lebed and DMS models do not 
break the paramagnetic limit H p , therefore the sharpen- 
ing of the anomaly in k(H) in dirty samples with lower 
H C 2 does not find natural explanation. Another appar- 
ent difficulty in the applicability of the LY theory is the 
existence of an additional phase boundary within the SC 
state, on which the increase of n starts in clean sam- 
ples. Yet, this point may be not so clear, since a cascade 
of phase transitions is predicted in Q1D case in DMS 
modeEj, some of which may survive in Q2D. From the 
above consideration, we see that FFLO theory gives a 
better description of our experiment, both in respect of 
the onset temperature and the tendency for the first- 
order transition in dirty samples. 

As it can be seen none of the features of our data con- 
tradicts formation of the FFLO state. For identification 
of this state more clear observation of the first order tran- 
sition in dirty samples may be usefulEil, however, final 
judgement can be made based on either phase-sensitive 
experimental, or direct observation of the structure of 
the order parameter in STM experiments. 

V. CONCLUSION 

In conclusion, the thermal conductivity in the organic 
superconductor A-(BETS)2GaCl4 in the magnetic field 
applied parallel to the Q2D plane shows a clear anomaly 
of the field dependence near H c2 , correlated with the 
shape of the H c2 \\ phase diagram at low temperatures. 
Both features disappear with the temperature increase 
and field inclination to the plane. This is consistent with 
the Pauli paramagnetic limiting in the low quality sam- 
ples, while with the FFLO state formation in the high 
quality samples. This observation may be of importance 
for explanation of the unusual field-induced superconduc- 
tivity in the closely related A-(BETS) 2 FeCl 4 saltElM 
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FIG. 1. Temperature dependence of the resistivity for 
A-(BETS) 2 GaCl 4 , samples #1 and #2. 
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TABLE I. 
A-(BETS) 2 GaCl 4 
marked with *. 


The properties of 6 single crystals 
Samples showing resistance jumps 


Sample No. 


# 1 


#2 


#3 


#4 


#5 


# 


batch 


A 


B 


A 


A 


B 


C 


p RT (fiilra) 


148 


190 


142 


170 


153 


14 


Prt/po 


150 


20 


23* 


3.5* 


1* 


13( 


Tc^mid (K) 


5.2 


5.6 


5.5 


5.1 


5.2 


5.: 


AT C (K) 


1.7 


2.2 


1.8 


1.6 


2.3 


1.! 


(-/ T )(T max) 
(«/ T )(T c ) 


1.35 


1.01 


1.38 


1.30 


1.06 


1.4 


ff c2 ||(0.3X) (T) 


12.7 


11.2 


12.9 


12.6 


11.1 


13. 


H*(0.3K) (T) 


9.7 


10.5 


9.7 


9.8 


10.7 


9.1 


Ak/k Hc2 (%) 


24 


6 


27 


32 


7 


34 



FIG. 2. Temperature dependence of k/T for 

A-(BETS) 2 GaCl 4 . Sample #1 at H=0 (curve 1-0) and in 
the normal state (1-1) {H= 6 T applied perpendicularly to 
the Q2D plane); sample #2 at H=Q (curve 2-0). 
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FIG. 3. Field dependence of k/T for sample 4f 
temperature of 0.3 K in a magnetic field perpendici 
plane (along b' direction, crosses), and parallel to 
along (triangles) and perpendicular (circles) to the 
direction (c axis). Double-ended arrow shows esti 
K e /T from Wiedemann-Franz ratio. 
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FIG. 4. Field dependence of k/T for samples #1 (a) and 
#2 (b). The data under the field parallel to the Q2D plane 
at several temperatures are shown by solid lines. The data 
for the perpendicular field at 0.3 K (shown with open circles) 
are shifted downward to avoid overlapping. 
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FIG. 5. Field dependence of k/T for high quality sa 
#1 and #6 and low quality samples #2 and #5. The 
were taken under the field parallel to the Q2D plane a1 
temperature 0.3 K. The curve for #6 is shifted upwar 
0.04 Wm^'K -2 and for #5 by the same amount down 
to avoid overlapping, the curves for #3 and #4 are simiL 
that for #1 and are not shown. 
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FIG. 6. Field dependence of p for 8 = (thin line) and 
of k under the fields inclined to the Q2D plane (thick lines, 
shifted upward with Q to avoid overlapping). T=0.3 K. (a) 
Sample #1, 6= 0°, 1°, 2 ° and 7° (from bottom to top), (b) 
Sample #2, 6= 0°, 1° and 5 °. 
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FIG. 7. Phase diagram for sample #2 c 
the midpoint of the resistive transition at fi 
down-triangles correspond to parallel and pe: 
orientations) and from the slope change in , 
(open circles). The inset shows the low-tem 
the phase diagram determined from k(H): 
circles are H c2 \\ for #1 and #2, respective! 
correspond to H* (see text) of #1. 
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FIG. 8. Dependence of fields H c2 \\ (circles) and H* (tri- 
angles) for #1 on the orientation of magnetic field within the 
conducting plane. The scale of the magnetic field is repre- 
sented in the lefthand side in units of tesla. The size of the 
symbol corresponds to error bar. The solid line is a guide for 
the eye. 
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